Spintronics-based mesoscopic heat engine 
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We consider a nanowire suspended between two spin-polarized leads and subject to a nonuniform 
magnetic field. We show that a temperature drop between the leads can significantly affect the 
nanowire dynamics. In particular, it is demonstrated that under certain conditions the stationary 
distribution of the mechanical subsystem has Boltzmann form with an effective temperature which 
is smaller than the temperature of the "cold" lead; this seems rather counter-intuitive. We also 
find that a change of the direction of the temperature gradient results in generation of mechanical 
vibrations rather than heating of the mechanical subsystem. 

PACS numbers: 73.63.-b, 73.23.Hk, 75.76.+j, 85.85.+j 



Nanomechanical resonators are devices which are be- 
ing employed not only to develop new technological ap- 
plications, such as ultra-sensitive sensors [US], but also 
to shed light on fundamental questions such as the tran- 
sition from the classical to the quantum mechanical de- 
scription of macroscopic objects [5]. Investigation of a 
system where a mechanical degree of freedom controls 
the properties of a mesoscopic junction between two bulk 
leads is an important line of research in nanomechan- 
ics [BHH]. In such structures, the mechanical part may be 
considered as a nanoengine whose operation is controlled 
by the states of two bulk thermodynamic reservoirs. It is 
well known that a macroscopic mechanical system may 
be driven into cyclic motion if it is coupled to reservoirs 
held at different temperatures. A Stirling engine, oper- 
ating by cyclic compression and expansion of air or other 
gas and placed between hot and cold spaces, is one ex- 
ample. 

By decreasing the size of the heat engine to the 
nanoscale level, quantum mesoscopic effects come into 
play and determine the behavior of both the working sub- 
stance and the mechanical subsystem. This opens new 
possibilities for the operation of a heat engine; for in- 
stance, reduction or even suppression of the mechanical 
fluctuations — effective cooling of the mechanical subsys- 
tem. Recently, it was shown that suppression, leading to 
ground state cooling, may be achieved if reservoirs, pre- 
sented by normal or superconducting metal leads, are 
held at different electrochemical potentials [9HT2]. It 
was also demonstrated that a temperature drop between 
reservoirs can also generate this effect if one assumes a 
very special "three particles" interaction inside a junc- 
tion [I3j[l4]. In this letter we investigate heating, pump- 
ing and cooling of a mechanical mode in a realistic nano- 
j unction where only a "two particle" interaction between 
mechanical degree of freedom and working subsystem ex- 
ists. We show that a temperature drop between linked 
leads can generate cooling or excitation, depending on its 
direction, of the mechanical subsystem. 

To be specific, we consider a carbon nanotube sus- 
pended between Lai-^Sr^MnOs ferromagnets of oppo- 
site polarizations (a structure recently realized experi- 
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FIG. 1: (a) A nanotube suspended between two spin- 
polarized leads and in proximity to a magnetic STM tip 
with magnetization M. The leads have opposite polariza- 
tions along the z-direction. (b) Nonuniform magnetic field, 
B, created by the magnetic tip. The nanotube (circles) de- 
flection, u, is in the ^-direction, (c) A doubly spin-degenerate 
electronic level exists in the nanotube at energy e - The ap- 
plied magnetic field splits this level into two levels a z =t, i 
separated by an energy A oc B z /jb- The leads are held at 
different temperatures Tl,r with a zero bias voltage. 

mentally [T5]) and subject to a nonuniform magnetic field 
B (cf. Fig. lb). The field can be generated, for example, 
by a magnetic STM tip in the form of a wedge which is in 
proximity to the nanotube (cf. Fig. la). If the nanotube 
is not too long (length L < 1/mi), it can be considered 
as a quantum dot coupled to the electrodes through tun- 
nel junctions. We assume that inside the nanotube there 
is only one doubly degenerate, with respect to spin, and 
spatially quantized electronic level which participates in 
electron exchange with the ferromagnetic leads. 

If the nanotube is straight and positioned below the 
STM sharp end, as shown in Fig. lb, the applied mag- 
netic field is directed towards the nanotube (z-direction) 
and Zeeman-splits the degenerate electronic level of the 
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nanotube. A two-level system (TLS) is thus formed in 
the latter with levels o z =t, 4- and energies e^, respec- 
tively. Deflection of the suspended part of the nanotube 
in the x-direction generates interlevel transitions (spin- 
flip) in the TLS. Consequently, the nonuniform magnetic 
field induces coupling between mechanical and electronic 
subsystems. In what follows, we will refer to this mech- 
anism of interaction between mechanical and electronic 
subsystem as spintronic coupling. 

We assume that the left and right leads are completely 
polarized along the z-direction (i.e., the density of states 

are ^^l) = ^ anc ^ V R(L) ^ Ur ( l ) ^ ^) and the energy 
difference between levels of the TLS is A = — > 
0. If there is no spintronic coupling and the intrinsic 
relaxation time r of the TLS is larger than the dwell- 
time of an electron in the nanotube, then the occupation 
number of the spin up (down) state inside nanotube is 
n T(l) - -M( e tU) _ e FL(R))/T L ( R) ) = fij(R) i where f F 
is the Fermi-distribution function, €fl(r) and ^l(h) are 
the Fermi energy and temperature of the left (right) lead, 
respectively. 

The spintronic coupling generates spin-flip transitions 
between the energy levels of the TLS. These transitions 
are inevitably accompanied by absorption or emission of 
mechanical quanta. The cooling (pumping) process of 
the mechanical subsystem is the result of an electron 
transition from the lower (upper) energy level to the up- 
per (lower) energy level of the TLS. In order to effec- 
tively cool (pump) the mechanical mode, it is required 
that <C n_i (n^ ^> n±). This can be realized by apply- 
ing a bias voltage or by subjecting the leads to different 
temperatures. Below, we consider cooling and pumping 
of the mechanical subsystem generated only by a tem- 
perature gradient between the leads. 

To perform a quantitative analysis of the system de- 
scribed above, we consider the Hamiltonian 

H = H nw + Hi + H t , (1) 
H nw = hujPb + e a h a + gu(d^di + (2) 

Hi = y^€L(fc)gfe.t.L a fc,t,£ + e -R(fc) a fe.l.fl a fc,+,-R> ( 3 ) 
k 

H t = ^Mt it>L ^ + ^4^fl^ + ft - c -» ( 4 ) 

k 

where CLk,a,L(R)(^l; a L(i?)) anc ^ ^<y(o\) being the annihila- 
tion (creation) operators for electrons in the left (right) 
leads and in the nanotube, respectively, and h a = a^d G . 
The first term in Eq. ([2| describes the nanotube mechan- 
ical degrees of freedom, which we restrict to the funda- 
mental flexural mode. This mode is described as a sim- 
ple harmonic oscillator with vibrational frequency uj and 
b(tf) being the annihilation (creation) operator for an ele- 
mentary excitation (vibron). The second term in Eq. ([2| 



describes the TLS with energy levels e a = cq ± A/2, 
where A is the Zeeman-splitting energy proportional to 
the z-component of the applied magnetic field B Z1 and 
eo is the zero-field energy, which is measured relative to 
the leads Fermi energy. The last term in Eq. Q de- 
scribes the spintronic coupling. It is proportional to the 
oscillator displacement, u = (fet + b)/y/2, and to the 
spin- flip operator. The coupling parameter g is equal 
to CfiBXod x B x (0) where jub is the Bohr magneton, xq is 
the zero-point vibrational amplitude, d x B x (0) is the field 
gradient along the x-direction, and C is a numerical fac- 
tor ~ A/L, where L is the nanowire length and A is the 
length of the STM wedge. This factor accounts for the 
electronic state inside the nanotube being extended over 
the whole length of the nanotube while the magnetic field 
is concentrated only in the region below the STM. The 
field gradient d x B x (0) induced by a Fe-based magnetic 
tip 7.5 nm thick with magnetization M = 1.75 • 10 6 A/m 
at a distance 7.5 nm is of the order of 35 mT-nm -1 . For 
this value of field gradient and for a nanotube with vi- 
brational frequency uj = 2tt • 100 MHz, mass 1 ag and 
A/L = 0.1, g is of the order of 2tt • 10 6 Hz. The term H x 
in Eq. (p3| describes the left and right leads. The term H t 
in Eq. H| describes the tunneling of electrons from the 
nanowire to the leads, and v.v., and are tunneling 

amplitudes [T6] . 

To analyze the performance of the system, we start 
from the Liouville-von Neumann equation for the total 
density operator g and then eliminate the lead electronic 
degrees of freedom [17]. If the temperatures of the leads 
are much greater than hw/k^^ the Fermi-distributions 
/l, r{z) ar e smooth functions within the energy interval 
huj. As a result one gets the following Lindblad master 
equation ([5| for the reduced density matrix p = Ttr+lq. 
The latter describes the mechanical degree of freedom 
and the electronic state of the TLS of the nanotube. 



d t p= --[H nw ,p] + C L [p) +C R \p], (5) 

where 

£>a[p] = r a [(l - f a )d aa pdl a + f a (dl a pd aa p) - 

-(1/2- f a ){n ak ,p}]. (6) 

Here, a = (L, R), cr L(i?) =t (I) and T a = 2n\t a \ 2 v a /h 
are the tunneling rates, v a is the density of states, and 
{A, B} denotes anticommutator. The collision integrals 
C a [p] describe the decoherence in the electronic subsys- 
tem induced by the bulk electronic reservoirs. 

If the resonant condition A = Huj is fulfilled and uj ^> 
g/h, T^, the rot at ing-wave- approximation can be used to 
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FIG. 2: Effective temperature T eff of the vibrational mode 
as function of the ratio G = T R /T L . For 6 > 1: T eff < 
mm{T L ,T R }. For G* < 9 < 1: T eff > max{T L ,T R }. For 
G < G* : no stationary distribution exits unless additional dis- 
sipation mechanisms are included. Inset: Vibron stationary 
distribution in logarithmic scale when the mechanical subsys- 
tem interacts only with ferromagnetic leads (curve 1), only 
with a bosonic bath at temperature Tb and coupling param- 
eter 7 (curve 2) and, both the fermionic and bosonic baths 
(curve 3). We use cj/2tt = 100 MHz, T L = 0.02 K, G = 10, 
T b = (T L + T R )/2, eo = k B T R /2, F L = V2g/H, F R = 1.62g/H 
and 7 = 0.002g/h. 



obtain the following rate equations 

P m (n) = ~g(P^(n)-Pi(n+l)), (7) 
P t (n) = -gP^(n+l) - T R f R P t (n) +T L f L P (n) 

-r L (i - f L )Pt(n) + r«(l - f R )P 2 (n), 
P 4 (n) = gP l n (n) +T R f R P (n) -T L f L P x {n) 

+r L (l - /L)P 2 (n) - T R {1 - f R )Pi(n), 
A(n) = -T L (l-f L )P 2 (n)-T R (l-f R )P 2 (n) 
+T L fLPdn)+T R f R P t (n), 
2Pi(n) = ~gn{P t (n - 1) - P+(n)) - (r L + r fl )^(n), 

where g = gV2/h. Here Pcr(^), ^2( n ) are the joint 

probabilities to find a vibrational mode in a Fock state 
with n vibronic quanta and one electron on the nanotube 
with spin a, empty nanotube and two electrons on the 
nanotube, respectively. Therefore, the total probability 
to find n vibronic quanta is P m (n) = Po(n) + P^(n) + 
Pl(n) + P2(n). The off-diagonal elements of the density 
matrix P l n (n) = lm(0|d t & (n_1) p(&t)^Jt| > / ( n - 1)! de- 
scribe the quantum entanglement between the electronic 
and mechanical subsystems generated by the correlation 
between spin-flip of electron and a change in the number 
of vibronic quanta. 

Equation ^ always has a stationary solution where 
P^(n) has a Boltzmann form 



P^(n) = Z _1 exp ( - Hun/k B T eff ), 



(8) 



FIG. 3: Vibron stationary distribution P m (n) when G < G*. 
We include a bosonic bath with temperature TJ, and coupling 
parameter 7. The stationary distribution exhibits a peak at 
rimax, which increases as 7 gets smaller. We use ribjh/g ^> 
1 (curve 1), ribjh/g = 0.021 (curve 2) and rib^h/g — 0.01 
(curve 3). Inset: The width of the distribution P m (n) scales 
inversely with T b . T b (A) = 1.5 K, T b (B) = 0.3 K, T b (C) = 
0.06 K and rib^h/g = 0.01 for all curves A, B and C. 



temperature is given by 



neff 



e-i + ^(i + e- 1 )!" 1 , (9) 



where Z = [l — exp(— fyjo /k^T e ^)\ 1 and the effective 



where 6 = Tr/Tl. This solution has physical mean- 
ing only if T e ff > 0. Negative effective temperature 
indicates that there is permanent energy pump into the 
mechanical subsystem. In this case in order to achieve a 
stationary regime one has to introduce additional exter- 
nal sources of dissipation. 

From Eq. ([9| one can see that for > and Tr > 
Tl, the vibrational mode is effectively cooled to a final 
temperature which is smaller than the temperature of 
the "cold" lead (T e *f < mm{T L ,T R } = T^) — cooling 
regime. From this equation also follows that in the inter- 
val 1 > 6 > 9* = (l-ficj/2e )/(l + fi^/2eo) the effective 
temperature is greater than the temperature of the "hot" 
lead (T e ff > max{Ti,,TR} = Tl) — heating regime. 

Note that when a temperature difference is used to re- 
duce mechanical fluctuations of the nanotube, the mini- 
mum effective temperature of the mechanical subsystem 
is always greater than T^/ n = 2T^/(1 + 2eo/huj), ap- 
proaching this value as — >• 00. As a result, the aver- 
age vibron number (n) = ^2 n TiP m (n) ~ exp(— co/feeTc) 
is finite but exponentially small. This is a consequence 
of the finite temperature of the "cold" lead which leads 
to infrequent pumping processes and does not allow to 
achieve absolute ground state cooling. 

From Eq. with T L and T R such that T e *f < 
(6 < 6*), it follows that d(n)/dt > and the aver- 
age number of vibrons increases with time. However, 
in a real physical situation the mechanical subsystem 
is also coupled to the phononic thermal baths of the 
leads. To account for the dissipation due to this coupling 
we add to the left side of Eq. ([5| a Lindblad operator 
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7>C 7 [19 . Here, the temperature of the phononic bath is 
Tb ~ (Tl + Tr)/2 and 7 = uo/Q is the coupling parame- 
ter characterized by the resonator quality factor Q, which 
can be very large Q ~ 10 5 at low temperatures [20 j. We 
properly modify the system of equations ([7]), and solve it 
numerically. The resulting shape of P^ (n) is depicted 
in Fig. [3j From this figure one can see that at small n 
the probabilities P^(n) increase with n. This is because 
for small n < T/g and 7 < g, the pumping rate (ex gn) 
is larger than the dissipation rate (ex 771). Then, the 
pumping rate eventually saturates to Y ~ minT^^) at 
n = n s ~T/g while the the dissipation rate governed by 
interaction with phonon reservoirs continue to increase 
linearly with n and finally overcomes the pumping rate. 
As a result, the distribution function reaches a maximum 
at n = Umax ~ r/7 and then exponentially decays for 
larger n > n max . 

Finally, we discuss conditions for ground state cooling 
of the mechanical subsystem; i.e., a stationary regime 
with a final vibron number (n) < 1. Thus far, we have 
considered the situation of completely polarized leads. 
However, in order to analyze the conditions for ground 
state cooling, it is necessary to estimate the effect of par- 
tial spin-polarization in the leads. To quantify the de- 
gree of polarizability of the right lead, we introduce a 
parameter t]l(r) — 1 — anc ^ ^ nc ^ that, for a 

symmetric case (r# = Tl), 1 — tjr <C 1, Tl <C Tr and 
eo ~ TrJzb ^> fouj, the average number of vibrons is 

(n) w (n) v=0 + (1 - f R )(l - m ) + 0((1 - m ) 2 ). (10) 

Hence, for a 90% spin-polarized right lead, the average 
vibron number is increased only by ~ 0.1. 

A necessary condition to achieve ground state cool- 
ing comes from Eq. ^ which suggests that in order 
to achieve the minimum effective temperature one needs 
Tr ^> Tl and eo ^> hu. The sufficient condition comes 
from requirement that the cooling rate k should be larger 
than the damping rate 7 due to a bosonic bath. The 
largest cooling rate for a given temperature gradient is 
achieved when > g/h > Tr(1 — Jr). Thus, the con- 
ditions for ground state cooling of the mechanical mode 
in the presence of a bosonic bath are Tr ^> Tl, eo ^> hu 
and Q ^> Jworib/gfR, where = [exp(/kj//cB^l>) — 1] _1 
(Tb is the temperature of the bosonic bath). From the 
last inequality, it follows that the best cooling regime is 
achieved when k^TR « eo- Our analysis shows that for 
a nanotube with frequency u = 2tt • 100 MHz (hw/k^ ~ 
6 mK) and realistic coupling parameter g « 2tt • 10 6 Hz 
(see above) and quality factor Q = 10 5 , the average 
vibronic number can be reduced to (n) = 0.44 for 
7^=200 mK and T L =20 mK. 

In conclusion we have studied heating, pumping and 
cooling of the mechanical vibrations of a nanotube sus- 
pended between two highly polarized magnetic leads. We 
have shown that spintronic coupling between the mechan- 
ical and electronic subsystems generated by a nonuniform 



magnetic field may result in suppression or generation 
of mechanical vibrations when the leads are held at dif- 
ferent temperatures. In particular, it was demonstrated 
that under certain conditions the stationary distribution 
of the mechanical subsystem has Boltzmann form with 
effective temperature which is smaller than the temper- 
ature of the "cold" lead. This counter-intuitive result 
is a consequence of the Fermionic nature of the ther- 
mal baths coupled to mechanical subsystem. Notice that 
coupling to bosonic baths results in an effective temper- 
ature of the vibrational mode equal to the mean value 
of the baths temperature. Also changing direction of the 
temperature gradient results in generation of mechanical 
vibrations rather than heating of the mechanical subsys- 
tem. Finally, for partial spin polarization in the leads, 
ground state cooling of the mechanical vibration can be 
achieved at realistic physical parameters if the leads have 
> 50% spin polarization. 
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